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Abstract. The roles of critical valence fluctuations of Ce and Yb are discussed as a key origin 
of several anomalies observed in Ce- and Yb-based heavy fermion systems. Recent development 
of the theory has revealed that a magnetic field is an efficient control parameter to induce 
the critical end point of the first-order valence transition. Metamagnetism and non-Fermi 
liquid behavior caused by this mechanism are discussed by comparing favorably with Celrlns, 
YbAgCu4, and YbIr2Zn2o. The interplay of the magnetic order and valence fluctuations offers 
a key concept for understanding Ce- and Yb-based systems. It is shown that suppression of the 
magnetic order by enhanced valence fluctuations gives rise to the coincidence of the magnetic- 
transition point and valence-crossover point at absolute zero as a function of pressure or magnetic 
field. The interplay is shown to resolve the outstanding puzzle in CeRhlns in a unified way. 
The broader applicability of this newly clarified mechanism is discussed by surveying promising 
materials such as YbAuCu4, /?-YbAlB4, and YbRh2Si2. 



1. Introduction 

In this article, we discuss the roles of critical valence fluctuations of Ce and Yb as a key origin of 
several anomalies observed in the Ce- and Yb-based heavy fermion systems. Recent development 
of the theory and experiments have revealed that the valence fluctuations appear ubiquitously 
by tuning the control parameters of the magnetic field and pressure, giving rise to broader 
consequences than previously recognized. First, we briefly survey accumulated experiments 
on the anomalies related to the valence transition and its fluctuations, and explain how the 
theoretical understanding has been achieved so far. Then, we discuss the recent development of 
theory and experimental results focusing on the prototypical materials. 

Valence transition is an isostructural phase transition with a valence of the materials element 
showing a discontinuous jump. A typical example is known as the 7-a transition in Ce metal [T], 
where the first-order valence transition (FOVT) occurs in the temperature-pressure (T, P) 
phase diagram while keeping the fee lattice structure, as shown in Fig. [IJa). The critical end 
point (CEP) of the FOVT is located at {T,P) =(600 K, 2 GPa). The valence of Ce changes 
discontinuously between Ce"'"^'''^ (7 phase) and Ce"*"^'^^ (a phase) at T = 300 K [2]. YbInCu4 
is also well known as a prototypical material for the isostructural FOVT [SI where as T 
decreases, the Yb valence changes discontinuously from Yb^^'^'^ to Yb^^'^^ at T = 42 K [5]. 
As diverging density fluctuations in the liquid-gas transition, valence fluctuations diverge at 
the CEP of the FOVT. When the temperature of the CEP is suppressed by controlling the 
materials parameters and enters the Fermi degeneracy regime, diverging valence fluctuations are 
considered to be coupled to the Fermi-surface instability. This multiple instability is considered 



to be a key mechanism for understanding anomalies observed in Ce- and Yb-based heavy- 
fermion systems. Such phenomena have been detected in CeCu2Ge2 [6] and CeCu2Si2 [7], where 
remarkable enhancement of the superconducting transition temperature appears at the pressure 
Pv larger than Pc corresponding to the antiferromagnetic (AF) quantum critical point (QCP) (see 
Fig. mb)), and in CeCu2(Sia;Gei_j;)2 where two separate domes of the superconducting phase 
have been found under pressure p|. In these materials, near P = P^, remarkable anomalies 
such as the T-linear resistivity and enhancement of the residual resistivity have been observed, 
which clearly exhibit distinct behavior of resistivity near P = P^ [HI HE]- Furthermore, these 
materials show that the coefficient A obtained by the fitting to the lowest-T regime of the 
resistivity decreases by 2-3 orders of magnitude when P increases across P = P^. 

Theoretically, it has been pointed out that enhanced Ce-valence fluctuations are the possible 
origin of the anomalies near P = P^ \n the CeCu2(Si/Ge)2 systems [S]. It has also been pointed 
out that the sudden drop of the coefficient A ai P = P^ can be understood as a sharp change of 
the valence of Ce on the basis of the Gutzwiller arguments in the periodic Anderson model |10j. 
It has been also shown that enhanced valence fluctuations cause the T-linear resistivity [7] 
and enhancement of residual resistivity [11]. The emergence of the superconducting dome near 
P ~ i-*v in Fig- has been explained by the valence-fiuctuation-mediated superconductivity 
shown by the slave-boson mean-field theory taking account of the Gaussian fiuctuations in the 
periodic Anderson model [10]. The density-matrix renormalization group (DMRG) calculation 
applied to the same model in one spatial dimension has also shown the enhancement of the 
superconducting correlation in the same parameter regime [12]. Pairing symmetry of density- 
fiuctuation (i.e., valence-fluctuation)-mediated superconductivity has also been analyzed on the 
basis of the phenomenological model [13j . 

Furthermore, recent NQR measurement in CeCu2Si2 has detected change of the Cu-NQR 
frequency near P = P^, suggesting that the Ce valence changes near P = P^ ^M- Detailed 
measurement of the T dependence of the specific heat and the upper critical field in CeCu2Si2 
under pressure concluded that the pairing symmetries of the superconducting phases near 
P = Pc and P = P^ are different, suggesting that a pairing mechanism different from the 
antiferromagnetic spin- fluctuation mediated one is realized near P^ |l^j . For a detailed summary 
of the theory and related experiments up to this stage, readers can refer to Ref. [16J and also 
Ref. [n]. 

Recent development of the theory of the quantum critical end point (QCEP) of the FOVT has 
revealed that valence fluctuations play a key role in the other Ce- and Yb-based heavy fermion 
systems more ubiquitously. We discuss the crucial roles of valence fluctuations in the following 
sections. In §2, we discuss the properties of the valence transition and valence crossover of Ce or 
Yb by comparing theoretical phase diagrams with Ce- and Yb-based systems. In §3, we discuss 
the fact that the QCEP of the FOVT can be induced rather easily by applying the magnetic 
field to the Ce- and Yb-based systems, which causes various anomalies such as metamagnetism 
and non-Fermi liquid behavior. We discuss that this mechanism gives a unified explanation for 
Celrlns, YbAgCu4, and YbIr2Zn2o. In §4, we discuss the interplay of magnetic order and valence 
fiuctuations, giving rise to coincidence of the magnetic transition and sharp valence crossover 
at T = K, as shown in Fig. [II|c). We show how the interplay resolves the outstanding puzzle 
in CeRhlns and also in YbAuCu4. The summary and outlook are given in §5. The broader 
applicability of the present mechanism is discussed by surveying promising materials such as 
/3-YbAlB4 and YbRh2Si2. 




Figure 1. (color online) (a) T-P phase diagram of Ce metal. The first-order valence transition 
line (bold line) terminates at a critical end point (filled circle). Schematic T-P phase diagrams 
for (b) Pc < Pv, and (c) Pc ^ Pv The dashed line represents valence-crossover temperature T* 
(see text). In (b) and (c), shaded regions represent superconducting phases. 

2. Model and phase diagrams 

2.1. Minimal model for Ce- and Yb-based heavy fermions exhibiting valence transition 

We consider the simplest minimal model, which describes an essential part of the physics in Ce- 
and Yb-based heavy fermion systems exhibiting valence transition, as follows: 

n = + Hi + Hi,y^ + Hu,^, (1) 

where He = J^ka ^kCko-Cko- represents the conduction band, H{ = ei J2ia ™L + ^ '^if'^^ 

the f level ef and onsite Coulomb repulsion U for f electrons, -ffhyb = ^Z^io- (/il-CiCT + cJo-Zio-) 

the hybridization V between f and conduction electrons, and Hfj^^ = Ufc X^i^i ^i'^f the Coulomb 
repulsion Ufc between f and conduction electrons. The Hu^^ term is a key ingredient for explaining 
the valence transition as well as the various anomalies caused by enhanced Ce- or Yb-valence 
fluctuations [TOl [121 [13 [H [IZ] : The T- linear resistivity and residual resistivity peak have also 
been shown theoretically by the model Eq. ([1]) [71 lllj. 

2.2. Ground state and finite-T phase diagrams 

Figure [21[a) shows the schematic ground-state phase diagram for paramagnetic states drawn on 
the basis of the DMRG calculation in one-spatial dimension (d = 1) |12j and the calculation 
by the dynamical- mean- field-theory (DMFT) in infinite-spatial dimension (d = oo) \W\ applied 
to the model Eq. ([T|). Here, the filling n = (nf -|- nc)/2 with na = J2aJ2iLi{0'lcr^ia) /N for a=f 
or c with N being the number of lattice sites is set slightly smaller than 1 (half filling), which 
describes typical Ce- and Yb-based heavy-fermion metals. 

The FOVT line (solid line) terminates at the QCEP (filled circle) and sharp valence crossover 
occurs (dashed line). When ef is deep enough, the Kondo state with nf = 1 is realized. As 
fff increases across the FOVT and valence-crossover line, the mixed valence (MV) state with 
rif < 1 is realized At the QCEP, the valence fluctuation Xv = —dui/def diverges and 

at the valence-crossover line, enhanced valence fluctuation appears [121 [18]. We note that the 
slave-boson mean-field theory ior d = 1 [12] , d = 2 [21] , and d = 3 [22l [18] also gives essentially 
the same phase diagram as Fig. [2]^a). This is because the valence transition and its fluctuation 




Figure 2. (color online) Schematic phase diagrams of (a) Ef-C/fc plane at T = K and (b) 
T-Ef-Ufc space for finite c-f hybridization V and large U in the periodic Anderson model Eq. ([1]) 
(see text). The first order valence transition line (bold line) terminates at the quantum critical 
end point (QCEP) (filled circle) and sharp valence crossover occurs on the thick dashed line. 
In (b), the FOVT surface (dark surface) and the valence-crossover surface (light surface) are 
separated by the critical-end line (thick line) touched at T = K, which is the QCEP. In (a) and 
(b), shaded regions at T = K represent superconducting phases. The thin dashed lines with 
arrow(s) in (a) represent routes for applying pressure to Ce metal (A) and Ce-based compounds 
(B), respectively. Corresponding T-P phase diagrams are shown by the cutout (A) and (B) in 
(b), which correspond to Fig. [TJa) and Fig. [T]^b) (or Fig. [l](c)), respectively (see text). 

are caused by the atomic Coulomb interaction Ufc in Eq. ([1]), which is ascribed to the local 
origin. Hence, the essential feature of the phase diagram for the valence transition does not 
depend on the spatial dimension. As shown by the shaded region in Fig. E^a), near the QCEP of 
the FOVT, the superconducting phase is shown to appear by the slave-boson mean field theory 
taking into account the Gaussian fluctuations [TO] , which is supported by the DMRG calculation 
in the same model Eq. ([1]) in d = 1 [12\ . 

Here, we remark that the valence transition is completely different from the localized-to- 
itinerant transition of f electrons. Everywhere in the phase diagram in Fig. [2][a), f electrons are 
itinerant with finite c-f hybridization (/j^^Cktr) 7^ 0. Hence, the number of f electrons is always 
included in the total volume of the Fermi surface, giving rise to the large Fermi surface (see 
also Fig. mb)). This is naturally understood in terms of the Landau-Luttinger sum rule. The 
existence of the QCEP in the ground-state phase diagram ensures that adiabatic continuation 
holds between the Kondo state and the MV state by detouring the QCEP of the first-order- 
transition line. 

Figure EJ^b) shows the schematic phase diagram in the T-ef-Ufc space. The FOVT surface 
(dark surface) bends to the side of the MV regime for at least large Ufc since larger entropy 
coming from f-spin degrees of freedom in the Kondo regime earns the larger free-energy gain [IBj . 
The FOVT line in Fig. E^a) is the line at the bottom of the FOVT surface. The critical-end 
line exhibiting diverging valence fluctuations, i.e., Xv = 00, separates the FOVT surface and the 
valence-crossover surface (light surface), which touches the T = plane at the QCEP. We note 
that even at the valence-crossover surface, valence fluctuations are enhanced [18], which cause 
remarkable anomalies in physical quantities in Ce- and Yb-based materials, as will be shown in 
53 and §4. 



2.3. Correspondence to Ce- and Yb-based materials 

As shown in Fig. [2)^a), the Kondo (MV) state with nf ~ 1 (nf < 1) is reahzed for deep (shallow) 
fff for a fixed C/fc. In Ce systems, the nf = 1 state corresponds to Ce"'"'^''' with a 4f ^-electron 
configuration per Ce site. In Yb systems, where Yb"''^''^ has 4f^^-electron configuration per Yb 
site, the hole picture is useful, so that the nf = 1 state corresponds to Yb'*'^'' with a 4f^-hole 
configuration per Yb site. 

When we apply pressure to the Ce (Yb) compounds, increases (decreases), since negative 
ions approach the tail of the wavefunction of 4f electrons (holes) at the Ce (Yb) site. The 
hybridization V and inter-orbital Coulomb interaction C/fc also increase. In the case of Ce metal, 
4f and 5d electrons are located at the same Ce site, which makes Ufc large [Ij. In Fig. [21^ a), the 
dashed line with an arrow (A) represents the route for applying the pressure to Ce metal. Hence, 
the cutout for the large C/fc including the T axis labeled by (A) in Fig. Wib) corresponds to the 
T-P phase diagram of Ce metal shown in Fig. [^a). Namely, the large critical-end temperature 
T ~ 600 K in Ce metal is naturally understood from a large Ufc due to its on-site origin in 
Fig.EJb). 

On the other hand, in the case of Ce (Yb) compounds with conduction electrons supplied 
from elements other than Ce (Yb), the Coulomb repulsion between the 4f electron at the Ce 
(Yb) site and the conduction electron Ufc is weaker in general because of its inter-site origin. 
Hence, applying pressure to the Ce compounds corresponds to the route represented by the 
dashed line with arrows (B) in Fig. [2]^a). The cutout for such a moderate Ufc labeled by (B) 
in Fig. El^b) is considered to correspond to the T-P phase diagrams of most of the Ce-based 
compounds. Although the magnetically-ordered phase is not shown in Figs. [2|^a) and EJ^b) for 
simplicity of explanation here, the cutout (B) corresponds to Fig. mb) (or Fig. mc)). Namely, in 
CeCu2(Si/Ge)2 systems, the AF phase is located in the Kondo regime with a certain interval to 
the valence-crossover pressure P^ in Fig. Hfb), where the valence-crossover temperature T*{P) 
covered by the superconducting dome corresponds to the cutout of Fig[2|^b). When we take 
into account the magnetic order in the phase diagram, depending on the strength of the c-f 
hybridization V (and also Ef and Ufc) in Eq ([TJ, the location of the magnetic-paramagnetic phase 
boundary changes. Hence, the relative position of Pc and P^ changes as shown in Figs. [UJb) 
and djc). The interplay of the magnetic order and valence crossover (or transition) is quite 
important in understanding the actual phase diagrams of Ce- and Yb-based systems, which will 
be discussed in §4. 

Here, we remind the readers of the fact that even at the 7-a transition in Ce metal, the 
magnitude of the valence jump is about 0.10 [2j as it is in the case of YbInCu4 [5], as remarked 
in §1. We note that the magnitude of the valence change at the valence crossover for moderate 
Ufc is smaller than that at the FOVT for large Ufc, which is expected to be in the order of 
0.01. Actually, such a tiny change of the Yb valence has been observed at the valence-crossover 
temperature T = ~ 40 K in YbAgCu4 [35|, which will be discussed in §3.3. 

3. Field-induced valence crossover 

3.1. Theoretical results 

As shown in §2, the valence transition is essentially ascribed to the charge degrees of freedom 
and it is nontrivial how the magnetic field affects the valence transition. Recently, theoretical 
studies have clarified that the QCEP of the FOVT as well as the sharp valence crossover is 
induced by applying a magnetic field [18]. The slave-boson mean-field theory applied to the 
model Eq. ([1]) in d = 3 with the Zeeman term —hJ^iiSf^ + Sf^) {h = gix-^H) has shown that 
the FOVT surface and the valence-crossover surface move as shown in Fig. [3)^ a) [22] I18j . On 
the locus with an arrow illustrated in Fig. [3)^ a), the magnetic susceptibility x = dm/dh with 
= T,i{fI^M - fljii + cJt'^it ~ diverges. Namely, metamagnetism occurs. The 

DMRG calculation applied to the same model in d = 1 has also shown that the field-induced 




Figure 3. (color online) (a) Magnetic-field dependence of FOVT surface and valence-crossover 
surface in the T-e^-Ufc phase diagram. The FOVT surface (dark surface) and the valence- 
crossover surface (light surface) are separated by the critical-end line (thick line) touched at 
T = K, which is the QCEP. The thick line with an arrow indicates the locus of the QCEP 
under magnetic field j22j. The black circle indicates a possible location of YbAgCu4 which 
touches the valence-crossover surface at about T = 40 K (black dashed line), (b) Ground-state 
phase diagram in ef-C/fc plane. The thick line with an arrow represents the locus of the QCEP 
under a magnetic field. The black square and circle indicate possible locations of Celrlns and 
YbAgCu4, respectively (see text and also Ref. [38j). The black star indicates a possible location 
of YbAuCu4. 

extension of the QCEP of the FOVT to the MV regime actually occurs, giving rise to the 
metamagnetisim [22]. This finding is quite important and hence we put special emphasis on this 
result: At the QCEP of the FOVT, not only the valence fluctuation, i.e., charge fluctuation, hut 
also the magnetic susceptibility diverges. 

Furthermore, this result indicates that even in the intermediate- valence materials, which does 
not show any valence transition at = 0, the QCEP of the FOVT is induced by applying the 
magnetic field H. An important point here is that as shown in Fig. [2jb), most of the Ce- and 
Yb-based compounds are located in the region for a moderate ?7fc, but not in the region for a 
large C/fc causing the FOVT. Hence, in the Ce- and Yb-based compounds, the valence-crossover 
temperature T*{H) with strong valence fiuctuations can emerge in the T-H phase diagrams. 
Here, the magnitude of the characteristic field to make T*{H) emerge depends on how close the 
location of the material is to the QCEP of the FOVT in the phase diagram. When the material 
is located close to the QCEP, metamagnetism as well as non-Fermi liquid behavior is expected 
to be observed by applying even a small magnetic field. In the following sections of §3.2, §3.3, 
and §3.4, we discuss the possible relevance of this mechanism to experimental observations. 

3.2. Metamagnesim and non-Fermi liquid behavior in Celrln^ 

Celrlns is a heavy-fermion metal, which shows a superconducting transition at T = 0.4 K at 
ambient pressure ^Sj- When pressure is applied, the superconducting transition temperature 
increases in spite of the fact that the In-NQR relaxation rate (TiT)~^ drastically decreases, 
suggesting the possibility that a pairing mechanism different from the spin-fiuctuation-mediated 
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Figure 4. (color online) T-K phase diagram in Celrlns [29j. For T < TpL, Fermi-liquid behavior 
appears in specific heat data [28]. 

one is relevant to Celrlns [23]. Interestingly, a metamagnetic anomaly was found to exist in the 
magnetization curve in Celrlns [25]. Several experimental groups have reported this anomaly 
occurring at the crossover (or first-order-transition) line in the T-H phase diagram, as shown 
in Fig. m [261 113 [2H1 [29]. As H increases to approach H^, non Fermi-liquid behavior becomes 
prominent: (i) Convex behavior appears in the T^'^ plot of the low-T resistivity suggesting 
that the T-linear resisivity appears near H = [28]. (ii) Residual resistivity has a peak at 
H = [29j . (iii) The specific-heat constant shows a logarithmic divergence C/T ~ — log T near 
H = |28j . The properties of (i) and (ii) are consistent with the theoretical results of critical 
valence fluctuations in Ref. [^ and Ref. [11], respectively, as explained in §1. The property of 
(iii) is also consistent with a theory recently developed for quantum valence criticality, which 
has shown that evaluation of the quasiparticle self-energy for one valence-fluctuation exchange 
process gives the low-T specific heat as C/T ~ — logT in a certain-T regime [30]. In Fig. 21 
the crossover line T*{H) touches T = K at H = ~ 28 T |29j . This fact and the above 
observations of (i)-(iii) suggest that Celrlns is located at the valence-crossover regime indicated 
by the black square in Fig. [3j^b). At H = 0, the black-square point is located at the Kondo 
regime. However, as H increases, the slope of the valence crossover line (dashed line) becomes 
steeper, which crosses the black-square point at ~ 28 T. This causes emergence of the valence- 
crossover line T*{H) for H > 28 T, as shown in Fig. HI Namely, Celrlns seems to be located at 
a distance of about H = Hy ~ 28 T from the valence-crossover line (dashed line in Fig. O^b)), 
close to the QCEP at H = 0. Hence, H = at T = K in Fig. H] seems to be in close proximity 
to the QCEP, which affects the low-T physics in Celrlns. 

Recent measurement of the de Haas-van Alphen (dHvA) effect has confirmed that the dHvA 
frequencies are not substantially altered at H = in Fig. |3] [2U] . This implies that the Fermi- 
surface topology remains essentially the same between the H < regime and the H > 
regime. This is also consistent with the metamagnetism caused by the field-induced QCEP of the 
FOVT: As explained in detail in Fig. [2]^a) in §2.2, the Fermi-surface volume basically does not 
change at the FOVT and valence crossover so long as a system remains in a paramagnetic state 
(Of course, in case of the FOVT (valence crossover), the lattice constant shows a discontinuous 
(continuous) change, which gives a change of the wave number itself). We also note that no 
evidence of the folded Brillouin Zone below and above has been obtained by the dHvA 
measurement in Celrlns p9]. This result indicates that the field-induced AF-ordered phase is 



unlikely to be realized at H^^ although the possibility of the magnetic-breakdown effect in the 
dHvA measurement should be carefully examined. 

We also note that Celrlns and CeCoIns have almost the same Fermi surfaces according to 
the band-structure calculations ^31j and the dHvA measurements ^32j. Both materials have a 
quite similar level scheme of the crystalline electric field (CEF) [33]. However, in CeCoIns, the 
crossover line (or first-order transition line) accompanied by the non-Fermi liquid behavior as 
shown in Fig. S] does not appear in the T-H phase diagram. We note here that by applying 
pressure to CeCoIns, residual resistivity and the coefficient A of the low-T resistivity drop 
drastically [33]. These observations suggest that the QCEP of the FOVT, or the sharp valence 
crossover point exists on a slightly-negative pressure side in the T-P phase diagram of CeCoIns. 
Actually, the experimental fact of the emergence of the T-linear resistivity and logarithmic 
divergence of the specific heat C/T ~ — logT at low temperatures and ambient pressure in 
CeCoIns does not contradict this point of view [30] . 

These results indicate that a viewpoint of the closeness to the QCEP of the FOVT is 
important for understanding the Cell5 systems in addition to the conventional view based on 
the competition between the Kondo effect and the RKKY interaction. Indeed, this viewpoint 
offers us a key to resolving the outstanding puzzle in CeRhlns as well, which will be discussed 
in §4.1. Since the superconducting correlation has been shown to develop near the QCEP of 
the FOVT theoretically (see Figs. [Sj^a) and[3];b)) [TOllIS], the proximity to the QCEP seems to 
control the occurrence of the unconventional superconductivity in Celrlns. 

To examine our theoretical proposal, it is desirable to measure the valence change of Ce at 
the crossover line T* (or first-order transition line) in Fig. 21 The X-ray adsorption spectra for 
direct observation of the Ce valence, the NQR measurement for the electric-field gradient, the 
ultrasonic measurement, and the X-ray diffraction and/or thermal expansion measurement for 
the lattice-constant change are highly desirable. 

3.3. Enhanced magnetic susceptibility and metamagnetism in YhAgCu4, 

As noted in §1, YbInCu4 shows the FOVT at T = = 42 K [3 S]. When In is replaced by 
other elements, the FOVT has not been observed. However, anomalous behavior which seems 
to be related to valence fiuctuations has been observed in YbAgCu4 [35]. In YbAgCu4, uniform 
magnetic susceptibility has a broad peak at T = 40 K. Below T = 40 K, the volume expansion 
occurs [36] simultaneously with a sharp valence crossover from Yb"*"^'^^ for T > 40 K to Yb"*"'^ '^^ 
for T < 40 K [35j, indicating that negative volume expansion occurs as T increases to approach 
T = 40 K. Namely, the uniform magnetic susceptibility is enhanced at the valence-crossover 
temperature T = = 40 K with strong valence fluctuations. This is consistent with our 
theoretical result discussed in §3.1 [18]. At the QCEP and the critical-end line in Fig. ^a), 
the valence fiuctuations diverge. At the same time, the magnetic susceptibility diverges [22j. 
YbAgCu4 seems to be located in the valence-crossover regime indicated by the black circle in 
Fig. [31^ a), whose distance from the valence-crossover surface at -ff = is about T = 40 K as shown 
by the black dashed line in Fig. [3l^a) [iSillTj . Hence, the broad peak of the magnetic susceptibility 
is considered to be caused by the valence fiuctuations developed at the valence-crossover surface. 
YbAgCu4 shows a metamagnetism in the low temperature limit around = 40 T |35] . This can 
also be naturally understood from the field-induced valence-crossover surface shown in Fig.[3]|^a). 
Namely, the location of YbAgCu4 indicated by the black circle in Fig. [Sj^b) is at a distance of 
about ~ 40 T from the valence-crossover line close to the QCEP at -ff = 0. For more details, 
the readers can refer to Refs. [17] and [18]. 

3.4-. Enhanced magnetic susceptibility and metamagnetism in YhIr2Zn2o 

Recently, behaviors similar to those of YbAgCu4 have been observed in YbIr2Zn2o [37]. The 
uniform magnetic susceptibility has a peak at T = = 7.4 K accompanied by the volume 



expansion for T < T*. Since the unit-cell volume with the Yb^^ '^ (4f'^ hole configuration) 
state is larger than that with the Yb^^ '' (4f ^ hole configuration) state, this observation suggests 
that the enhancement of the magnetic susceptibility is caused by the enhanced Yb-valence 
fluctuations associated with the valence change of Yb. This viewpoint is consistent with the fact 
that the thermal expansion coefficient is negative for T <T* j37) . 

Recently, it has been found that YbIr2Zn2o shows a metamagnetism at ~ = 10 T [37]. 
Near the metamagnetic field H^, enhancement of the cyclotron mass of electrons has been 
observed by the dHvA measurement. The specific-heat constant C/T and the coefficient of 
the resistivity at low temperatures also indicate the mass enhancement of electrons at H = H^. 
Furthermore, residual resistivity has a peak at H = H^. The temperature region where the 
low-T resistivity shows the dependence becomes narrowest in the vicinity of H = H^, 
which suggests a tendency of the T- linear resistivity near H = H^. All these observations 
are naturally explained by the field-induced valence crossover discussed in §3.1. Although the 
mass enhancement and the residual-resistivity peak suggest a clear signature of the Yb-valence 
fiuctuations developed when H approaches H^, the experimental fact that the resistivity has the 

dependence at the lowest temperatures even at H = implies that YbIr2Zn2o at ambient 
pressure is located at the valence-crossover regime. Namely, YbIr2Zn2o seems to be located at 
a similar position to that of YbAgCu4, indicated by the black circle in Fig. [3]^b). The distance 
from the valence-crossover line close to the QCEP at H = seems to be about ~ 10 T in 
YbIr2Zn2o, since metamagnetic field is about ~ 10 T [38]. Indeed, a steep volume shrinkage 
has been observed at ~ [37] , which indicates that the Yb valence sharply increases at 
H = H^. The dHvA measurement has concluded that the Fermi surfaces do not change at 
H = from the fact that the dHvA frequencies do not change across [37J. This is also 
consistent with the mechanism of the field-induced valence crossover, since the Fermi-surface 
volume is essentially unchanged at the valence crossover so long as the system remains in the 
paramagnetic state as discussed in §2.2 and also in §3.2. 

When pressure is applied to YbIr2Zn2o, the metamagnetic field can be tuned to approach 
if = T |39j. This seems to be consistent with the location of YbIr2Zn2o mentioned above, since 
applying pressure makes C/fc large because of the reduction of the distance between wavefunctions 
of 4f and conduction electrons and also makes < small (i.e., |ef| large) in the hole picture. 
Here, we should note that the panel of the cutout shown in Fig. EJ^b) is drawn for the Ce 
compounds. In the case of Yb compounds, the panel tilts with a certain angle with the opposite 
sign to the ef axis. The observation of the quantum criticality at the QCEP of the FOVT 
by tuning the control parameters of the pressure and/or magnetic field is an interesting future 
issue [40J. Systematic measurements in YbX2Zn2o (X=Ir, Co, and Rh) are highly desirable for 
a unified understanding of the effect of the critical valence fluctuations. 

4. Interplay of magnetic order and valence fluctuations 

In the preceding sections, we have discussed the nature of the QCEP of the FOVT in a 
paramagnetic phase and how it is controlled by applying the pressure and/or magnetic fleld. 
In Ce- and Yb-based heavy fermion systems, however, the magnetic order can also occur. 
The interplay of the magnetic order and Ce- or Yb-valence fluctuations offers a key concept 
for understanding these systems. Below we focus on CeRhlns as a prototypical material to 
demonstrate how such an interplay plays a crucial role in understanding the puzzling behaviors. 

4.1. CeRhln^ 

CeRhlns is a heavy-fermion metal j41j which undergoes an AF transition at Tn = 3.8 K with 
the ordered vector Q = (1/2, 1/2,0.297) at ambient pressure [IS]. When pressure is applied, the 
AF-ordered phase changes to the paramagnetic and superconducting phase at P = Pc ~ 2 GPa, 
as shown schematically in Fig. [5] [431 H5] H6l ITT] HS] H9] . This material has attracted much 
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Figure 5. (color online) Schematic T-P phase diagram of CeRhlns [431 144] , Shaded region 
represents superconducting phase. 

attention since accumulated experiments offer outstanding puzzle, whose significant features are 
summarized as follows: (i) The Sommerfeld constant 7e ~ 50 mJ/(molK-^) in the AF state at 
P = is about 10-times larger than 7e = 5.7 mJ/(molK^) in LaRhlns |41l 146) . (ii) The Fermi 
surfaces similar to those in LaRhlns in the AF phase for P < P^ ~ 2.35 GPa change to the 
Fermi surfaces similar to those in CeCoIns in the paramagnetic phase for P > Pc by the dHvA 
measurement performed under the magnetic field ~ 15 T |46j. (iii) The cyclotron mass of 
electrons shows an enhancement toward P = Pc- The cyclotron mass of the /32 branch, whose 
Fermi surface has a cylindrical shape, changes from drriQ at P = to 60mo at P < Pc and the 

signal is not detected for P > Pc, probably because of a too heavy mass, ~ lOOmo [46j. (iv) The 
T- linear resistivity emerges prominently near P = Pc [431 [17] . (v) Residual resistivity has 
a peak at P ~ Pc [lU HU [JT]. (vi) The superconducting phase exists in a wide pressure region 
around P = Pc [ISl [H [35] . 

Because of the dHvA measurement noted above (ii), one might succumb to the temptation 
to believe a scenario in which the localized to itinerant transition of f electrons happens at 
P = Pc [501 I51j . However, this scenario encounters a serious difficulty in explaining the 
experimental fact (i) above: The AF state with the 10-times mass enhancement at P = 
strongly suggests that the heavy quasiparticles contribute to the formation of the AF state, 
indicating the existence of the c-f hybridization even in the AF phase for P < Pc. Hence, the 
theoretical explanation for resolving a series o/ puzzles outlined above (i)-(vi) in a natural and 
unified way has been desired. 

First of all, let us point out that the transport anomalies above (iv) and (v) and the 
robust superconducting phase (vi) are quite similar to the observations in the CeCu2(Si/Ge)2 
systems O [Zl IS] introduced in §1. Emergence of the T- linear resistivity prominent at P = Pv 
in Fig. [It^b), at which the residual resistivity has a peak, covered with the superconducting 
phase in the wide pressure region is the common feature in CeCu2Ge2 [6], CeCu2Si2 [Z\, and 
CeCu2(Gea;Sii_x)2 [S]. Since the T-linear resistivity and residual resistivity peak appear at the 
AF-paramagnetic boundary in CeRhlns, Pc ~ Py seems to be realized in Fig. [SJ As will be 
shown below, this is actually the case realized in the model Eq. ([T]) for realistic parameters for 




Figure 6. (color online) (a) Ground state phase diagram in e^-Ufc plane of the model Eq. ([T]) for 
V = 0.2 and U = oo at n = 0.9 (see text) [21j . The FOVT (solid line with triangles) terminates 
at QCEP (filled circle). Valence crossover with enhanced valence fluctuations occurs at the 
dashed line with open circles. The solid line with filled squares represents the antiferromagnetic 
(AF)-paramagnetic (PM) boundary. Occupied (solid lines) and empty (gray lines) bands in the 
periodic Anderson model at n = (rif + nc)/2 = 0.9 are shown in the paramagnetic phase (b) 
and in the AF-ordered phase (c) . In (c) , the gray dashed line indicates the energy band of the 
conduction band, ek at ric = 0.8. The black dashed line is a guide for the eyes indicating that the 
Fermi surface in AF-ordered phase with finite c-f hybridization coincides with the small Fermi 
surface where f electrons are completely localized. 

CeRhlns, which naturally resolves the above puzzles (i)-(vi) |21] . 

Figure El^a) shows the ground-state phase diagram of the model Eq. ([1]) without the magnetic 
field on the square lattice obtained by the slave-boson mean-field theory [52], which enables us 
to treat the AF order and valence transition or fluctuations on equal footing. Here, we set small 
hybridization V = 0.2 and the conduction band as = — 2t(cos kx + cos ky) at the total filling 
n = 0.9 to mimic the /32 branch of CeRhlns [l6|. The FOVT line (solid line with triangles) 
terminates at the QCEP (filled circle). Valence crossover with enhanced valence fiuctuations 
occurs on the dashed line with open circles. These are the results within the paramagnetic 
states. When the AF order with the ordered vector Q = (vr,7r) is taken into account, we find 
that the AF-paramagnetic boundary represented by the solid line with filled squares emerges at 
almost the same place as the FOVT line and the valence-crossover line. Since £f can be regarded 
as the pressure as explained in §2.3, this result implies that Pq ~ Pv is realized. 

The precise explanation of this mechanism is illustrated in Figs. El^a) and[7[b). As shown 
in Fig. [T[a), when the Neel temperature close to K, i.e., ?n(-P) ~ K, meets the valence- 
crossover temperature T*{P), enhanced valence fluctuations suppress the AF order. Then, the 
AF order is suddenly cut around P = P^, as shown in Fig. [7|^b). Actually, the slave-boson 
mean-fleld calculation shows that the AF order exhibits the flrst-order transition at the solid 
line with filled squares in Fig. EJa) ^21j. Namely, the coincidence of the AF transition and the 
valence crossover or transition is caused by the enhanced valence fiuctuations or the FOVT. 
Hence, the Pc ~ Pv illustrated in Fig. W^c), or more precisely in Fig. [7^b), is shown to be 
realized for a small hybridization case in the model Eq. ([1]), which is a realistic parameter for 




Figure 7. (color online) (a) Schematic T-P phase diagram for < P^- (b) Enhanced valence 
fluctuations at Pv suppress the AF order, giving rise to the AF transition of the first order at 
P = Py shown by vertical thick line [2T] . 

CeRhlns [l6j. We should remark that for larger hybridization V, the AF-paramagnetic boundary 
appears at smaller ef than the FOVT and valence-crossover line in Fig. [6)^a), which corresponds 
to Fig. mb). Hence, depending on the strength of the hybridization V, the relative position of 
Pc and Pv changes as shown in Figs, mb) anddJ^c), and Fig. El^a). 

We have also shown that the Fermi surface calculated for V = 0.2 in the AF phase at 
n = (nf + nc)/2 = 0.9 is nearly the same as the small Fermi surface without c-f hybridization 
where electrons for nf = 1 are located at the localized f level and the conduction band is filled 
up with all extra electrons for Uc = 0.8 |21] . This naturally explains the dHvA measurement 
(ii) above. In the paramagnetic phase in Fig. E](a), the large Fermi surface, which counts the 
f-electron number in the total Fermi volume, i.e., rif + ric, is realized, as shown in Fig. El^b). In 
the AF phase, the lower and upper hybridized bands are folded, as shown in Fig. E[c). Since the 
lower folded hybridized band is completely filled, the Fermi surface becomes the same as that of 
the conduction electrons at the filling ric- Hence, we stress that the small Fermi surface appears 
in the AF phase with finite c-f hybridization, {fl.^c\^„) 7^ 0. 

The mass enhancement observed by the dHvA measurement (iii) above is also quantitatively 
reproduced by the same model Eq. ([T]) under the magnetic field if = 15 T |21] . Namely, the mass 
enhancement of the two-dimensional-like Fermi surface of the j32 branch from 6mo at P = to 
60mo at P < Pc [46] is well reproduced by the present mechanism, which shows that the total 

density of states at the Fermi level increases near the AF-paramagnetic boundary in Fig. EJ^a). 
This is because when the pressure is applied to the AF state, i.e., when increases, the gap 
between the original lower hybridized band and the folded band increases as in Fig. [6)^c). Then, 
the f-electron-dominant flat part of the folded lower hybridized band approaches the Fermi level, 
giving rise to the increase in the density of states. In the paramagnetic phase, as ef approaches 
the AF-paramagnetic boundary in Fig. [6]^a), i.e., as ef decreases, nf increases to approach 1, 
which is in the so called Kondo regime. Hence, the density of states at the Fermi level increases. 
Thus, the reason why the mass enhancement occurs toward the AF-paramagnetic boundary 
P = Pc in CeRhIn5 is naturally explained. 

This result is quite consistent with a recent experiment on the eff'ective mass of electrons [43j : 
Knebel et al applied a magnetic fleld to CeRhlns, whose magnitude is the same as if = 15 T 
used in the dHvA measurement |46j. After the disappearance of the superconductivity, the 



dependence of the resistivity appears at lowest temperatures. Knebel et al found that the ^/A 
with A being the coefficient of the term scales nicely with the cyclotron mass of the P2 branch 
under pressure for < P < Pc, i.e., y/A/m* =const. scaling holds [l3]. This scaling indicates 
that the mass enhancement toward P = P^ is essentially caused by the effect of the energy band 
dispersion of quasiparticles [53]. We also note that about 10-times mass enhancement at P = is 
also reproduced by the same model Eq. ([1]) [21], which is quite consistent with the experimental 
fact (i) above. Formation of the heavy quasiparticles via the c-f hybridization naturally explains 
the mass enhancement even inside of the AF phase. 

As shown in Fig. [6|^a), Pc ~ Pv is a consequence of the suppression of the AF order by 
enhanced valence fluctuations for moderate C/fc, e.g., C/fc ~ 0.5, which is a realistic parameter 
for CeRhlns, giving rise to the T- linear resistivity and the residual resistivity peak at the AF- 
paramagnetic boundary. Hence, the experimental facts (iv) and (v) are naturally explained. 
Since the superconducting phase was shown to be realized around P^ [THl US], as shown in 
Fig. mb), the above (vi) is also consistently explained (see Fig. [5] and Fig. [Tl^b)). Hence, all of 
the above experimental observations (i)-(vi) are naturally explained in a unified way. 

We note that local correlation effects of 4f electrons give rise to the large ordered moment of 
the AF order: Only about 10 % reduction of the full ordered moment for the Kramers doublet 
of the lowest CEF level is realized even near the AF-paramagnetic boundary in Fig. [6l^a) by the 
slave-boson mean-field theory [21], which is consistent with the observation in CeRhlns [48j. Our 
result shows that strong onsite Coulomb repulsion of 4f electrons yields the local-moment-like 
behavior as if f electrons are localized, in spite of the fact that f electrons are always itinerant, 
(/ko-'^kcr) / always holds everywhere of the phase diagram shown in Fig. E^a), which is 
nothing but a manifestation of the localized-itinerant duality [5l] ■ 

As reported in Ref. [55], the experimental comparison between CeRhlns and 
CeCo(Ini_2;Cda;)5 shows that CeRhlns at P > Pc corresponds to CeCoIns at P = [16|. This 
is quite consistent with our argument in §3.2: Because Pc ~ P^ is realized in CeRhlns, if the 
above correspondence holds, that implies that on the slightly-negative pressure side in the T-P 
phase diagram of CeCoIns, the QCEP or sharp Ce-valence crossover point, i.e., P = Pv, exists, 
which is the same conclusion stated in §3.2. 

The interplay of the magnetic order and enhanced valence fluctuations plays an important 
role not only in CeRhlns but also in other materials. One such example is YbAuCu4, which will 
be discussed in the next section. 

4.2. YhAuCui 

YbAuCu4 is a sister compound of YbInCu4, which exhibits the FOVT as noted in §1. We 
note that the phase diagram of YbAuCu4 shown in Fig. [5] is essentially the same as Fig. (Tj^c) 
or Fig. [Tl^b), if P is replaced by H [56]. The ^^Cu-NQR measurement has detected that there 
exists a characteristic temperature T*(H) in the T-H phase diagram, at which the Cu-NQR 
frequency changes sharply, strongly suggesting that the Yb- valence crossover occurs ^7\. The 
Yb- valence crossover temperature T*{H) becomes = K at the field H^, which seems to 
coincide with the field at which the Neel temperature goes to Tn = K [56^ i57j (see Fig. [8|) . 
Namely, ~ occurs in the T-H phase diagram. As discussed in §3.1, the valence crossover 
temperature T*{H) is induced by applying the magnetic field. If YbAuCu4 is located at the 
moderate Uic regime indicated by the black star in Fig. ^h), the appearance of T*{H) in Fig. [8] 
is naturally explained |17j . Namely, YbAuCu4 seems to be in the Kondo regime at = and 
by applying H = ~ 1.3 T, the valence crossover T*{H) is induced, as shown in Fig. [8l This 
indicates that YbAuCu4 is close to the QCEP at -ff = with a distance of about ~ 1.3 T as 
indicated by the black star in Fig. [3j[b). As explained in §4.1, when T^(H) close to K meets 
T*{H), a suppression of the magnetic order by enhanced valence fiuctuations occurs, which gives 
rise to the coincidence He ~ H^- In YbAuCu4, the cusp-like anomaly of the residual resistivity 
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Figure 8. (color online) T-H phase diagram of YbAuCu4 [S^. The Neel temperature Tn 
is suppressed and the Yb valence crossover temperature T* (thick dashed line) is induced for 
H > ~ 1.3 T. For T < TpL (black dashed line), the resistivity shows the Fermi liquid 
behavior p ^ T^. 



was observed at H = He ~ [57]. The mass enhancement of electrons estimated from the A 
coefficient in the resistivity was also observed near H = ^ |57j . These observations are 
naturally explained by the enhanced Yb- valence fluctuations |17] . 



5. Summary and outlook 

In this article, the roles of critical valence fluctuations in Ce- and Yb-based heavy fermion metals 
are discussed. Recent development of theory and experiment has revealed that critical Ce- and 
Yb- valence fluctuations play a key role in several anomalies in this family of materials. Magnetic 
field is shown to be a useful control parameter to induce the quantum critical end point (QCEP) 
of the first-order valence transition (FOVT). As a prototypical material, metamagnetism and 
non-Fermi liquid behavior in Celrlns, YbAgCu4, and YbIr2Zn2o are shown to be naturally 
explained by the mechanism of the field- induced valence crossover. The interplay of the magnetic 
order and Ce- or Yb-valence fluctuations is a key mechanism for understanding anomalous 
behavior in Ce- and Yb-based heavy fermion systems. It is shown that the interplay resolves 
the outstanding puzzle in CeRhlns. The origin of the transport anomalies and drastic change 
of the Fermi surfaces accompanied by the mass enhancement of electrons is naturally explained 
by the interplay in a unified way. A special emphasis is put on the fact that a small Fermi 
surface generally appears by the folding of the hybridized band in the AF phase in the periodic 
Anderson model. The unified explanation for Celnlns and CeRhIn5 strongly suggests that the 
viewpoint of the closeness to the QCEP of the FOVT is indispensable in understanding the 
Cells systems. As another prototypical material where the interplay occurs, the T-H phase 
diagram of YbAuCu4 is discussed. The proximity to the QCEP of the FOVT is a key concept 
for understanding Ce- and Yb-based heavy fermion systems as an underling mechanism. 

For outlook, recently synthesized materials, YbIr2Zn2o, YbCo2Zn2o, and YbRh2Zn2o are 
promising materials, which are expected to show the properties caused by the Yb-valence 
fluctuations by tuning the control parameters of the magnetic field and pressure. Recently, in 
/3-YbAlB4, where the superconductivity was first discovered among Yb-based compounds [61j, 



evidence of strong Yb-valence fluctuations has been reported [62]. Further study to clarify 
the role of the valence fluctuations in the unconventional criticality as well as the origin of 
the superconductivity is an interesting future problem. We point out that the T-H phase 
diagram in YbRh2Si2 [58], where unconventional criticality has attracted much attention in 
correlated electron systems [59], is closely similar to that in YbAuCu4 (see Fig. [8]). Hence, there 
exists a possibility that the critical Yb-valence fluctuation plays a key role in the origin of the 
unconventional criticality in YbRh2Si2. Thus, it is quite important to examine the possibility 
of whether the Yb-valence change takes place at the crossover temperature T*{H) whose origin 
has not yet been clarifled experimentally in the T-H phase diagram. The experimental fact that 
YbRh2Si2 shows a huge mass enhancement even at H = such as 7e ~ 1.7 J/(molK^) deep 
inside of the AF phase [60j indicates that heavy quasi particles form the AF-ordered state as 
in CeRhlns. The Co-NQR measurement in Yb(Rhi_a;Co2;)2Si2 and/or the direct observation 
of the Yb valence by X-ray adsorption measurement at the crossover temperature T*{H) are 
highly desirable. As noted in §2.3, the magnitude of the valence change is expected to be of the 
order of 0.01 in most of the Ce and Yb compounds. Hence, high accuracy measurement which 
can detect the tiny change of the valence is a challenging future problem, which will open new 
avenues of study in this field. 
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